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ABSTRACT 

We present a new database of absorption and emission-line measurements based on the entire 
spectral atlas from the Sloan Digital Sky Survey (SDSS) 7*'' data release of galaxies within a redshift 
of 0.2. Our work makes use of the publicly available penalized pixel-fitting (pPXF) and gas and 
absorption line fitting (gandalf ) IDL codes, aiming to improve the existing measurements for stellar 
kinematics, the strength of various absorption-line features, and the flux and width of the emissions 
from different species of ionised gas. Our fit to the stellar continuum uses both standard stellar 
population models and empirical templates obtained by combining a large number of stellar spectra 
in order to fit a subsample of high-quality SDSS spectra for quiescent galaxies. Furthermore, our 
fit to the nebular spectrum includes an exhaustive list of both recombination and forbidden lines. 
Foreground Galactic extinction is implicitly treated in our models, whereas reddening in the SDSS 
galaxies is included in the form of a simple dust screen component affecting the entire spectrum that 
is accompanied by a second reddening component affecting only the ionised gas emission. In order to 
check for systematic departures from the rather standard set of assumptions that enters our models, 
we provide a quality assessment for our fit to the SDSS spectra in our sample, for both the stellar 
continuum and the nebular emissions and across different wavelength regions. This quality assessment 
also allows the identification of objects with either problematic data or peculiar features. We hope to 
foster the discovery potential of our database ; therefore, our spectral fit is available to the community. 
For example, based on the quality assessment around the Ha and [Nil] A6584 lines, approximately 1% 
of the SDSS spectra classified as "galaxies" by the SDSS pipeline do in fact require additional broad 
lines to be matched, even though they do not show a strong continuum from an active nucleus, as do 
the SDSS objects classified as "quasars". Finally, we provide new spectral templates for galaxies of 
different Hubble types, obtained by combining the results of our spectral fit for a subsample of 452 
morphologically selected objects. 

Subject headings: catalogs — galaxies: abundances — galaxies: stellar content — galaxies: statistics 



1. INTRODUCTION 

One of the most powerful ways to unveil the forma- 
tion and evolution of galaxies is to study their spectra. 
Since William WoUaston first observed sunlight pass- 
ing through a prism in 1802 to reveal a rich spectrum, 
spectral lines have been used to understand the physical 
processes governing the evolution of stars and galaxies. 
Galaxy spectra feature stellar absorption and nebular 
emission lines, both of which have been studied exten- 
sively over the past several decades. 

The prominent absorption features seen in galaxy spec- 
tra reflect the average surface gravities, effective temper- 
atures, metal abundances and kinematics of their stel- 
lar components and thus hold a preserved record of a 
galaxy's formation and evolution. In order to measure 
the strength of various stellar absorption lines so as to 
compare them across different objects and facilitate the 
prediction of theoretical models we require a common 
system of integrated light spectroscopy, the most promi- 
nent of which is the Lick/IDS system (Burstein et al. 
1984; Gorgas et al. 1993; Worthey et al. 1994; Worthey 
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& Ottaviani 1997). This widely- used system spans the 
optical wavelength regime from CN at -4150A to TiO, 
at ^6200A and defines a total of 25 standard absorp- 
tion line indices. Measurements of these absorption line 
strengths, combined with stellar population models have 
been effectively used to study galaxy properties, such as 
: luminosity-weighted ages, metallicities and abundance 
ratios of early-type galaxies (Davies, Sadler & Peletier 
1993; Fisher, Franx & lUingworth 1995, 1996; Jorgensen, 
I. 1999; Kuntschner 2000; Kuntschner et al. 2006). 

Nebular emission lines, on the other hand, probe the 
physical state of the ionized gas in galaxies and thus 
can be used to trace the nuclear activity of, for exam- 
ple, a central supermassive black hole or their instanta- 
neous rate of star formation (Osterbrock & Pogge 1985). 
Emission-line ratios have been turned into powerful di- 
agnostic tools, not just for individual galaxies, but for 
massive spectroscopic surveys. Baldwin, Phillips & Ter- 
levich (1981) first proposed the use of diagnostic dia- 
grams, which have subsequently been used and refined by 
numerous research groups (Veilleux & Osterbrock 1987; 
Kewley et al. 2001; Kauffmann et al. 2003; Schawinski et 
al. 2007). 

Throughout the past two decades the Sloan Digital Sky 
Survey (SDSS) has established the largest and most ho- 
mogeneous database of both photometric and spectro- 
scopic data for galaxies. Using the Apache Point Ob- 
servatory's 2.5 m telescope, the SDSS has obtained 9400 
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deg^ of imaging in the u, g, r, i and z bands ( Fukugita 
et al. 1996; Gunn ct al. 1998, 2006) on the AB magnitude 
system (Smith et al. 2002), considering both astrometric 
(Pier et al. 2003) and photometric (Ivezic et al. 2004) 
calibrations. Follow-up spectroscopy (Eisenstein et al. 
2001; Strauss et al. 2002; Richards et al. 2002; Stoughton 
et al. 2002; Blanton et al. 2003) has furthermore yielded 
nearly one million objects spectroscopically classified as 
"galaxies" (See Tab. 1 and Fig. 1 of Abazajian et al. 
2009). The seventh, and latest, data release of the SDSS 
(SDSS DR7, Abazajian et al. 2009) is notably improved, 
in terms of both photometry and spectroscopy, and is to 
date the largest sized database. 

The SDSS spectroscopic database has been useful in 
numerous investigations; however, the SDSS pipeline 
measurements still suflFer from a few crucial shortcom- 
ings. The most serious is that the pipeline vahies for the 
absorption-line strength are not corrected for the impact 
of nebular emissions. Not accounting for emission infills 
leads to an underestimation of the true strength of var- 
ious absorption-line features and, in objects with strong 
emission lines - such as star-bursting systems or active 
nuclei - a more refined treatment of the nebular com- 
ponent may be needed to secure an unbiased extraction 
of the stellar kinematic. Furthermore, when measuring 
the ionised-gas emission it is desirable to include red- 
dening by dust and impose a sensible a priori on the 
relative strength of recombination lines, such as those 
from the Balmer series. The MPA-JHU group have ad- 
dressed some of these issues (e.g., Tremonti et al. 2004), 
refining the absorption-line strength measurements for 
SDSS galaxies and leading to many useful studies (e.g., 
Kauffmann et al. 2004; Cid et al.' 2005; Kewley et al. 
2006; Erb et al. 2006; Faber et al. 2007). Yet, in the case 
of both these databases, it is not clear how systematic 
problems may have affected the extraction of the phys- 
ical parameters listed in these catalogues, and it is not 
always possible to verify the quality of these measure- 
ments themselves. In this paper we attempt to further 
improve both the absorption and emission-line spectro- 
scopic measurements of nearby SDSS galaxies, as well 
as providing an assessment of the quality of our fit to 
both the stellar and nebular spectral components of our 
sample galaxies. This helps to isolate objects with prob- 
lematic data or systems that cannot be described by our 
set of physical assumptions and which may thus con- 
tain additional structures. For instance, in this paper 
we will show how such a quality assessment will reveal a 
previously unidentified population of objects with broad 
Balmer lines. Similar inspections may help finding the 
presence of other additional features such as the typical 
"bump" at 4640A associated to the presence of Wolf- 
Rayet stars (e.g., Brinchmann, Kunth & Durret 2008) or 
a systematic mismatch in the NaD region around 5900A 
due to the presence of neutral material entrained in galac- 
tic outflows (e.g., Sato et al. 2009, Chen et al. 2010). In 
fact, it is with the hope of fostering the potential discov- 
ery of such peculiar systems that we enable the commu- 
nity to inspect each of our fits with the SDSS galaxies 
considered in this study. 

This work, unlike the other databases, adopt both the- 
oretical stellar population synthesis model (Bruzual & 
Chariot 2003) and empirical stellar templates (Sanchez- 



Blazquez et al. 2006) to improve accuracy of the fit. Be- 
sides, we treat nebular emission which affects all over the 
optical wavelengths including [Ni] AA5198,5200 doublet 
that is not considered by MPA-JHU group. Dramatic 
effect of [Ni] AA5198,5200 doublet subtraction on Mgb 
absorption line strength is described in § 5.3. Specifically, 
our sample includes the entire SDSS DR7 spectroscopic 
database of objects that the SDSS pipeline identified as 
"galaxies" , which we restricted by imposing a redshift cut 
below z < 0.2. This yielded a total of 661,400 galaxies. 

This work is organized as follows. In § 2, we describe 
our method for deriving the absorption and emission 
line measurements, which includes extracting the cen- 
tral stellar and gas kinematics. In § 3 we explain how 
we assessed the quality of our fit in both the stellar con- 
tinuum and the nebular emission, and explain how to 
identify and deal with objects with broad Balmer lines. 
The structure of our new database is then described in 
§ 4, with § 5 containing the details of how our measure- 
ments for the stellar kinematics, the emission-line fluxes 
and widths, and the strength of absorption-line features 
compare with the same quantities from the SDSS DR7 
pipeline or the MPA-JHU DR7 catalogue. As an exam- 
ple of the improvements made by using our procedure, at 
the end of § 5 we discuss how accounting for the presence 
of the weak [Ni] AA5198, 5200 emission can impact esti- 
mates of the element abundances in galaxies. Finally, we 
present a number of stacked spectra for different Hub- 
ble types and discuss how they can be useful to galaxy 
studies in § 6 and give a summary of our work in § 7. 

2. METHOD 

In order to measure the stellar and gaseous kinematics 
within the central regions of our sample galaxies and as- 
sess the depth of the most important stellar absorption- 
line features, as well as the strength of the nebular emis- 
sion observed in the SDSS spectra, we adopted a strat- 
egy very similar to that used during the SAURON integral- 
fleld spectroscopic survey of nearby early-type galaxies 
(Emsellem et al. 2004; Sarzi et al. 2006; Kuntschner et 
al. 2006), which dealt with more than 34,000 individual 
spectra. 

Our analysis of the SDSS spectra consisted of three 
separate steps. First, we extracted the stellar kinemat- 
ics by directly matching the spectra with a set of stellar 
templates, while excluding the spectral regions poten- 
tially affected by nebular emission. We then measured 
the gaseous kinematics and assessed the strength of each 
emission-line by fitting the entire spectrum, simultane- 
ously using the stellar templates that were broadened 
and shifted by the previously derived stellar kinematics 
and a set of Gaussian templates representing the nebular 
emission. Finally, we subtracted, from the SDSS spectra, 
our best model for the nebular emission and measured 
the strength of various stellar absorption lines from the 
observed spectra following standard line-strength index 
dcifinitions. An early application and description of this 
fitting strategy can be found in Schawinski et al. (2007) 
and Cardamone et al. (2009). 

2.1. Stellar Kinematics 

To measure the stellar kinematics from the SDSS spec- 
tra in our catalogue, we used Cappellari & Emsellem's 
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TABLE 1 

lONIZED-GAS EMISSION LINES 



No 


Species Wavelength (A) 


relative strength to (line) 


1 


[On] 


3726.03 




2 


fOiil 


3728.73 




3 


[Neiii] 3868.69 




4 


Nciii 3967.40 




5 


H5 


3889.05 


037 CRrv) 


6 


He 


3970.07 


056 CHn) 


7 


H5 


4101.73 


091 (Un) 


8 


H7 


4340.46 


164 fHn-^ 


9 


[Olll 


4363.15 




10 


[Hen 


4685.74 




11 


[Ariv 


4711.30 




12 


[Ariv 


4740.10 




13 


H/3 


4861.32 


0.350 (Ha) 


14 


[Oiii 


4958.83 


350 (fOlIll A5007~) 


15 


[Oiii 


5006.77 




16 


[Nl] 


5197.90 




17 


[Ni] 


5200.39 




18 


[Hcl 


5875.60 




19 


[Oi] 


6300.20 




20 


[Oi 


6363.67 


0.333 ([Or] A6300) 


21 


[Nil 


6547.96 


0.340 ([Nil] A6584) 


22 


Ha 


6562.80 




23 


[Nil] 


6583.34 




24 


[Sii] 


6716.31 




25 


[Sii] 


6730.68 




26 


[Arm] 7135.67 





(2004) pixel-fitting method, Icnown as pPXF, and param- 
eterized the line-of-sight velocity distribution (LOSVD) 
by means of a simple Gaussian, thus obtaining only the 
galaxy redshift z and central stellar velocity dispersion 
fT*. During the pPXF fit, we used the Bruzual & Char- 
lot (2003) stellar population models used by Trcmonti et 
al. (2004), which we combined with a number of empiri- 
cal stellar templates based on the IVIILES stellar library 
(Sanchez-Blazquez et al. 2006). The use of the JVIILES li- 
brary helps with matching absorption-line features dom- 
inated by alpha-enhanced elements (see §2.4). Before 
extracting the central LOSVD in our sample galaxies, 
we also needed to define the spectral regions used during 
the pPXF fit. This means setting the wavelength range of 
the fit, excluding the spectral regions potentially affected 
by nebular emission, masking the regions where the sky- 
lines at 5577A, 6300A, and 6363A(in the rest frame) may 
have left substantial residuals after their sul)traction, and 
finally excluding the NaDAA5890, 5896 absorption lines 
that were not properly matched by the stellar templates 
owing to the impact of interstellar absorption. The wave- 
length range of the pPXF fit was set by the spectral cover- 
age of our templates, which is restricted to 68OOA owing 
to the presence of the IVIILES templates, and the redshift 
estimate provided by the SDSS for each object, which 
also helped with the placement of a 1,200 km s^^-wide 
mask (for emission-lines with a velocity dispersion of 200 
km s~^) centered around the expected position of each of 
the lines listed in Tab. 1. During the pPXF fit, we used a 
4*^-order additive polynomial correction for the spectral 
slope of our templates. 



2.2. Emission Line Measurements 



Once the stellar kinematics had been measured, we 
lifted the emission-line mask and proceeded to simul- 
taneously match the stellar continuum and the nebu- 
lar emission in our sample galaxies using Sarzi et al.'s 
(2006) gandalf code. During this matching process, 
gandalf combines the stellar templates with a number 
of emission-line templates represented either by single 
Gaussians, in the case of single emission-line species, or 
by sets of Gaussians meant to describe either doublets, 
such as [O iii]AA4959, 5007, or multiplets, such as the en- 
tire Balmer series. The intrinsic relative strength of the 
lines, prior to any reddening due to dust in the target 
galaxies, was set by either atomic physics, in the case of 
doublets, or by assuming a gas temperature, in the case 
of the Balmer lines. During the gandalf fit the position 
and the intrinsic width (prior to convolution with the 
line-spread fimction of the SDSS spectra) of the Gaussian 
templates were solved non-linearly by means of a stan- 
dard Levenberg-Marquardt minimization^, while at each 
iteration of this process the relative strengths of both the 
stellar and emission-line templates were solved linearly 
and by imposing non-negative weights. We imposed the 
same kinematics on all of the forbidden lines and the 
recombination lines. Only the gas kinematics was de- 
rived at this stage. During the entire gandalf fit the 
stellar templates were shifted and broadened according 
to the previously derived stellar LOSVD. Reddening by 
dust was included in the geindalf fit and was solved non- 
linearly along with the gas kinematics. We included two 
reddening components in our models: the first which af- 
fects the entire spectrum, representing dust diffusion ev- 
erywhere throughout the target galaxy; and the second, 
which affects only the nebular emission and can account 
for dust more localized around the emission-line regions. 
Finally, the NaD and skyline regions were excluded dur- 
ing the gandalf fit and, prior to that we de-reddened 
the SDSS spectra to correct for foreground Galactic ex- 
tinction (using maps from Schlegel et al. 1998) and to 
measure the intrinsic estimates for dust reddening in our 
sample galaxies. A complete list of the emission lines 
that were fit, including their intrinsic relative strengths, 
is given in Tab. 1. 

2.3. Absorption Line Measurements 

The last step of our analysis consisted of measuring 
the strength of various stellar absorption lines through- 
out the wavelength range of the SDSS spectra. This was 
done after subtracting our best gandalf model for the 
nebular emission and using the standard line-strength 
definition for the indices listed in Tab. 2. To perform the 
index measurements we used an IDL routine kindly pro- 
vided by H. Kuntschner. Prior to the index calculation, 
this routine brought the spectra back to the rest-franic 
and is capable of degrading their spectral resolution to 
match that of the Lick/IDS system, to which most of the 
indices listed in Tab. 2 belong. Whether our line-strength 
measurements were based on the original emission-line 
cleaned spectra or on their degraded counterparts (in 
our catalogue we provide both values), we corrected our 
index values for the impact of stellar kinematic broaden- 

For such non-Uncar fit both pPXF and gandalf use 
the MPFIT IDL routine provided by Craig Markwardt at 
http: //cow. physics. wise. edu/"craigm/idl/ 
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TABLE 2 
Absorption line indices 



No 


Name 


Index Bandpass (A) 


Pseudocontinua(A) 


Reference 


01 


H5a 


4083.500-4122.250 


4041.600-4079.750, 


4128.500-4161.000 


Lick'' 


02 


H(5f 


4091.000-4112.250 


4057.250-4088.500, 


4114.750-4137.250 


Lick 


03 


CNi 


4142.125-4177.125 


4080.125-4117.625, 


4244.125-4284.125 


Lick 


04 


CN2 


4142.125-4177.125 


4083.875-4096.375, 


4244.125-4284.125 


Lick 


05 


Ca4227 


4222.250-4234.750 


4211.000-4219.750, 


4241.000-4251.000 


Lick 


06 


G4300 


4281.375-4316.375 


4266.375-4282.625, 


4318.875-4335.125 


Lick 


07 


H7A 


4319.750-4363.500 


4283.500-4319.750, 


4367.250-4419.750 


Lick 


08 


H7f 


4331.250-4352.250 


4283.500-4319.750, 


4354.750-4384.750 


Lick 


09 


Fe4383 


4369.125-4420.375 


4359.125-4370.375, 


4442.875-4455.375 


Lick 


10 


Ca4455 


4452,125-4474.625 


4445.875-4454.625, 


4477.125-4492.125 


Lick 


11 


Fe4531 


4514.250-4559.250 


4504.250-4514.250, 


4560.500-4579.250 


Lick 


12 


C4668 


4634.000-4720.250 


4611.500-4630.250, 


4742.750-4756.500 


Lick 


13 


H/3 


4847.875-4876.625 


4827.875-4847.875, 


4876.625-4891.625 


Lick 


14 


H/3g 


4851.320-4871.320 


4S1 ^ nnn 4^/1^ nnn 


4ss;n nnn 4Q'^n nnn 


Gonzalez'^ 


15 


Fc4930 


4903.000-4945.500 


^on /iQO? nnn 




Gonzalez 


16 


OIIIi 


4948.920-4978.920 


ziRS^ nnn aol'x^ nnn 
^ooo.uuu-^yoo.uuu. 


f^n'^n nnn ^n7n nnn 


Gonzalez 


17 


OIII2 


4996.850-5016.850 


4885.000-4935.000, 


5030.000-5070.000 


Gonzalez 


18 


Fe5015 


4977.750-5054.000 


4946.500-4977.750, 


5054.000-5065.250 


Lick 


19 


Mgi 


5069.125-5134.125 


4895.125-4957.625, 


5301.125-5366.125 


Lick 


20 


Mg2 


5154.125-5196.625 


4895.125-4957.625, 


5301.125-5366.125 


Lick 


21 


Mgb 


5160.125-5192.625 


5142.625-5161.375, 


5191.375-5206.375 


Lick 


22 


Fe5270 


5245.650-5285.650 


5233.150-5248.150, 


5285.650-5318.150 


Lick 


23 


Fe5270s 


5256.500-5278.500 


5233.000-5250.000, 


5285.500-5308.000 


SAURON<= 


24 


Fe5335 


5312.125-5352.125 


5304.625-5315.875, 


5353.375-5363.375 


Lick 


25 


Fe5406 


5387.500-5415.000 


5376.250-5387.500, 


5415.000-5425.000 


Lick 


26 


Fe5709 


5696.625-5720.375 


5672.875-5696.625, 


5722.875-5736.625 


Lick 


27 


Fe5782 


5776.625-5796.625 


5765.375-5775.375, 


5797.875-5811.625 


Lick 


28 


NaD 


5876.875-5909.375 


5860.625-5875.625, 


5922.125-5948.125 


Lick 


29 


TiOi 


5936.625-5994.125 


5816.625-5849.125, 


6038.625-6103.625 


Lick 


30 


Ti02 


6189.625-6272.125 


6066.625-6141.625, 


6372.625-6415.125 


Lick 



Wortiicy ct al. 1994 

Gonzalez, J. J. 1993, PhD thesis, Univ. California, Santa Cruz 1993 
Kuntschner et al 2006 



ing. To achieve this effect, we observed how the values of 
the indices changed when we measured them on the opti- 
mal combination of the stellar templates returned by the 
gandalf fit {LS optimal) and on the best stellar model 
obtained during the same process (LSmodei), which is 
the optimal template convolved by the LOSVD obtained 
during the pPXF fit and further adjusted for dust redden- 
ing by gandalf. Given the index measurement on the 
emission-cleaned SDSS spectra (LS), the corresponding 
index value corrected for kinematic broadening can be 
derived as follows 

r C — T Q ^^optimal n 

Prior to the linci-strciigth measurements we subtracted 
all of the emission fitted by gandalf, even that from weak 
lines that would be too weak to be formally deemed de- 
tected. Although this may produce an over-subtraction 
of the nebular flux and result in somewhat biased values 
of the absorption- line indices, we prefer this to introduc- 
ing a sudden step in the index infill correction by not 
subtracting the emission when crossing a certain detec- 
tion threshold. Moreover, when the skylines fall within 
the pseudocontinuum or index passbands of one of our 
indices, which can happen for different indices at dif- 
ferent redshifts, we provided the LSopUmai value that is 
measured on the stellar optimal template. The residuals 
of an imperfect subtraction of such skylines may lead to 
meaningless values for the indices if these were measured 
directly on the SDSS data. 



2.4. Optimization of Template Library 

As briefly mentioned at the beginning of this section, 
the templates based on the stellar population models 
from Bruzual & Chariot (2003) cannot fully describe the 
stellar spectrum of our sample galaxies. This is partic- 
ularly true in cases of the most massive or the closest 
objects in our sample, where the SDSS spectra probe 
overall or central stellar populations that display an over- 
abundance of alpha-elements, which are not currently 
accounted for by the population models. Typical spec- 
tral regions that are not well matched in these cases are 
those around the CN 4200A band or the Mgb index. If 
stellar population models cannot fully match such alpha- 
enhanced galactic spectra, the combination of a large 
number of the stellar spectra on which the models are 
based can. Such stellar spectra cannot be used individu- 
ally together with the population models, as this would 
make the pPXF and gandalf fit excessively cumbersome. 
Furthermore, the resulting description of the stellar con- 
tinuum would not necessarily be physically motivated, 
which in turn could lead to biased emission-line results 
and a wrong infill correction for the absorption-line in- 
dices. On the other hand, following Sarzi et al. (2010) we 
could select a number of SDSS spectra of objects totally 
devoid of nebular emission to construct, with the stellar 
library, a more limited number of semi-empirical alpha- 
enhanced templates. More specifically, using an initial 
subsample of our catalogue, we drew a number of passive 
galaxies covering a range of values for the stellar velocity 
dispersion, to trace the galaxy mass, and for the ratio 
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of the Mgb and Fc5015 absorption-lino indices, to trace 
the degree of alpha-enhancement (sec, e.g., Thomas et al. 
2005). Using pPXF we then matched each of these spec- 
tra with the entire MILES stellar library from Sanchez- 
Blazquez et al. (2006), without masking the emission-line 
regions and while employing only multiplicative polyno- 
mials to adjust the shape of the stellar templates. The 
quality of such a fit is generally excellent, thanks to the 
large number of templates, to the extent that the re- 
sulting optimal combinations of the stellar spectra used 
in each of these fits could, in practice, be regarded as 
a galactic spectrum deconvolved from the line-of-sight 
kinematical broadening. Such optimal templates were 
then themselves combined into 12 empirical templates 
and designed to cover, as evenly as possible, the distri- 
bution of the fitted galaxies in the (T*-Mgb/Fe5015 space. 
Finally, these templates were convolved and re-binned in 
order to match the spectral resolution and sampling in 
the Bruzual & Chariot (2003) model templates. 

2.5. Errors 

Both the pPXF and gandalf returned formal errors for 
the stellar and gaseous kinematics and the flux values of 
the fitted emission lines. Prom MonteCarlo simulations 
based on synthetic spectra created from our best gandalf 
fit to a subsample of galaxies in our catalogue and the 
formal uncertainties of their SDSS spectra, we verified 
that the formal uncertainties were consistent with the 
scatter of the parameter estimations obtained from fit- 
ting the synthetic spectra. Errors on the line-strength 
indices for each single object in our database would be 
more cumbersome to compute. In fact, although the im- 
pact of the statistic uncertainties on the flux densities 
of the SDSS spectra could be easily described in terms 
of formal errors in the line indices, the knock-on effect 
of the uncertainties on the stellar kinematics and - most 
importantly - on the gandalf emission-line fit would re- 
quire a MonteCarlo simulation in order to be properly 
estimated. Since this is not viable for nearly 700,000 ob- 
jects, we decided instead to attach to any given line-index 
measurement in the objects of our catalogue the typical 
uncertainty that is measured in spectra showing a sim- 
ilar quality, a similar stellar kinematic broadening and 
strength of the emission falling in the pseudo-continuum 
of index passband of the considered index. To derive such 
typical errors, we selected a large number of SDSS spec- 
tra covering a wide range in quality (quantified by the 
average value of the S /N ratio between the values of the 
flux density and the statistical SDSS errors on them), of 
values of the stellar velocity dispersion cr*, and - within 
similar S/N and ct* values - of strength for the emission 
possibly affecting the line indices (e.g., [Om] A5007 for 
Fe5015). The strength of the emission lines was quan- 
tified by the ratio A/rN of the best fltting amplitude A 
of the considered line and the level of scatter rN in the 
residuals of the gandalf fit. Based on our best gandalf 
fit for these objects and the formal SDSS errors on the 
flux densities, we created a number of synthetic spectra 
with a slightly different kinematic broadening and neb- 
ular emission, as allowed by the formal uncertainties on 
the pPXF and gandalf parameters, and re-ran our en- 
tire analysis on them to obtain new values for the line- 
indices. By re-matching each of these synthetic spectra 
and subsequently measuring the line-strength index af- 



ter subtracting our best model for the nebular emission, 
we obtained a range of absorption-line strength values 
and the corresponding Icr-errors. The scatter of such 
synthetic line-index values around the one measured on 
the original spectra provided us with the typical errors 
on the line-strength measured in our catalogue. Fig. 1 
shows how their values vary as a function of S/N, cr^, and 
- when applicable, the A/N of the contaminating lines. 
Furthermore, despite the presence of nebular emission, 
the error budget in the absorption-line indices is chiefly 
dominated by the quality of the spectra themselves. The 
formal uncertainties are also provided for the values of 
both reddening components included in our gandalf flt. 
The typical errors for absorption-line indices are listed in 
Tab. 2.5. 

2.6. Final Remarks 

We conclude this section with a couple of remarks. 
First, before the strength of the absorption features were 
assessed at the end of our procedure, the first two steps 
of our analysis were repeated in order to optimize our 
gandalf fit. In particular, once the position and the 
width of the lines were measured the first time, we re- 
peated the pPXF fit by placing a better emission-line mask 
on only the region where some emission was found. This 
makes us more confident in the stellar kinematics, in 
particular for galaxies with little or no emission or in 
star-forming objects where the Balmer absorption lines 
are the predominant feature in the stellar spectrum and 
cannot be properly matched unless the mask around the 
Balmer emission lines - generally narrower than the ab- 
sorption lines - is properly adjusted. A better stellar 
kinematic leads to a better gandalf fit. 

A second point to keep in mind is that the second com- 
ponent for the dust extinction used in the gandalf fit - 
that affecting only the nebular spectrum - should only be 
considered reliable in the presence of Balmer lines, and 
thus of an observed Balmer decrement. If no Balmer lines 
are detected, or if only Ha is. then the de- reddened values 
for the fluxes of all the lines fitted by gandalf will also 
be very uncertain, and one should only use the observed 
fluxes. Finally, we stress that even in the presence of a 
Balmer decrement in a good quality spectrum, the val- 
ues of the diffuse component for the reddening - approx- 
imated as a uniform screen affecting the entire galaxy 
spectrum - may to some extent be degenerate with the 
particular gandalf choice for the stellar templates. De- 
termining the star-formation history and the exact val- 
ues of the reddening in our sample galaxies is beyond the 
scope of this catalogue, however, and here we are inter- 
ested mostly in providing the best physically motivated 
and accurate description of the stellar and nebular spec- 
trum, in order to provide reliable values for the stellar 
and gas kinematics, for the emission-line fluxes and for 
the infill-corrected absorption-line indices. 

3. FITTING QUALITY ASSESSMENT 

Assessing the quality of our fits to the SDSS spectra is 
a key to the accuracy of our measurements for the stellar 
and gaseous kinematics, as well as for the strengths of 
both the emission and absorption lines. An inadequate 
model for the observed spectrum, or artificial features in 
the data themselves is likely to introduce biases to the 
measured parameters that are not accounted for by our 
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TABLE 3 

Typical errors for absorption-line indices 



S/N''~ 10 S/N~ 20 S/N~ 30 

Name A/N'^<4 10<A/N<25 30<A/N<50 A/N<4 10<A/N<25 30<A/N<50 A/N<4 10<A/N<25 30<A/N<50 



HdA 


1, 


.434±0, 


,405 


1, 


,046±0, 


,212 





,859±0, 


,139 


0, 


,462±0 


,193 





,574±0, 


,100 


0, 


,439±0, 


,017 


0, 


,323±0, 


,088 


0, 


,332±0 


,018 





,375±0, 


,036 


H<5f 


0, 


.960±0, 


,294 


0, 


,700±0, 


158 


0, 


,557±0, 


,087 


0, 


,400±0, 


,066 


0, 


,379±0, 


,051 


0, 


,292±0, 


,040 


0, 


,247±0, 


,055 


0, 


,235±0, 


,018 


0, 


,214±0, 


,023 


CNi 


0, 


.038±0, 


,010 


0, 


,031±0, 


,007 


0, 


,028±0, 


,005 


0, 


,013±0, 


,001 


0, 


,018±0, 


,003 


0, 


,013±0, 


,001 


0, 


,009±0, 


,003 


0, 


,010±0, 


,001 


0, 


,010±0, 


,001 


CN2 


0, 


.051±0, 


,013 


0, 


,035±0, 


,007 


0, 


,030±0, 


,005 


0, 


,024±0, 


,006 


0, 


,019±0, 


,003 


0, 


,015±0, 


,001 


0, 


,014±0, 


,003 


0, 


,010±0, 


,000 


0, 


,012±0, 


,001 


Ca4227 


0, 


.540±0, 


,149 


0, 


,442±0, 


,091 


0, 


,382±0, 


,042 


0, 


,259±0, 


,049 


0, 


,258±0, 


,046 


0, 


,193±0, 


,012 


0, 


153±0, 


,022 


0, 


,168±0, 


,030 


0, 


,152±0, 


,016 


G4300 


1, 


.180±0, 


,309 


0, 


,905±0, 


159 


0, 


,832±0, 


,125 


0, 


,507±0, 


,073 


0, 


,494±0, 


,067 


0, 


,414±0, 


,065 


0, 


,283±0, 


,027 


0, 


,305±0, 


,028 


0, 


,297±0, 


,030 


H7A 


1, 


.345±0, 


,395 


1, 


,019±0, 


,221 


0, 


,899±0, 


,116 


0, 


,568±0, 


,082 


0, 


,559±0, 


,086 


0, 


,468±0, 


,067 


0, 


,352±0, 


,045 


0, 


,338±0, 


,040 


0, 


,435±0, 


,051 


H7f 


0, 


.866±0, 


,283 


0, 


,592±0, 


114 


0, 


,538±0, 


,087 


0, 


,356±0, 


,045 


0, 


,330±0, 


,045 


0, 


,264±0, 


,032 


0, 


,214±0, 


,031 


0, 


,184±0, 


,009 


0, 


,225±0, 


,020 


Fc4383 


1, 


.454±0, 


,324 


1, 


,242±0, 


,220 


1, 


,160±0, 


,129 


0, 


,599±0, 


,060 


0, 


,621±0, 


,107 


0, 


,571±0, 


,068 


0, 


,417±0, 


,044 


0, 


,387±0, 


,014 


0, 


,444±0, 


,046 


Ca4455 





.655±0, 


,183 


0, 


,559±0, 


,088 


0, 


,523±0, 


,066 


0, 


,315±0, 


,035 


0, 


,312±0, 


,037 


0, 


,257±0, 


,024 


0, 


181±0, 


,018 


0, 


,205±0 


,026 





,211±0, 


,028 


Fc4531 


1 


.048±0, 


,218 


1, 


,020±0, 


176 


0, 


,968±0, 


,103 


0, 


,493±0 


,064 





,511±0, 


,081 


0, 


,426±0 


,042 


0, 


,316±0, 


,025 


0, 


,313±0 


,003 





,328±0, 


,056 


C4668 


1, 


.912±0, 


,295 


1, 


,745±0, 


,356 


1, 


.742±0, 


,234 


0, 


,848±0, 


.068 


0, 


.881±0, 


,118 


0, 


,758±0, 


.054 


0, 


,493±0, 


,048 


0, 


,564±0, 


.047 


0, 


.578±0, 


,022 


H/3 


0, 


.726±0. 


,133 


0, 


,646±0. 


,266 


0, 


.623±0, 


,068 


0, 


,352±0, 


.027 


0, 


.299±0. 


,031 


0, 


,283±0, 


.029 


0. 


,190±0, 


,009 


0, 


,208±0, 


.004 


0, 


.232±0. 


,026 


H/3g 


0, 


.523±0. 


,110 


0, 


,460±0. 


,180 


0, 


.462±0, 


,049 


0, 


,242±0, 


.027 


0, 


.238±0, 


,037 


0, 


,213±0, 


.021 


0. 


,148±0, 


,009 


0, 


,140±0, 


.004 


0, 


.161±0. 


,012 


Fe4930 


0, 


.912±0. 


,384 


0, 


,885±0. 


,237 





.908±0, 


,105 


0, 


,539±0, 


.100 


0, 


.534±0. 


,231 


0, 


,401±0, 


.049 


0. 


,274±0, 


,025 


0, 


,315±0, 


.013 


0, 


.298±0. 


,012 


OIIIi 


0, 


.693±0. 


,147 


0, 


,743±0. 


,123 





.649±0, 


,127 


0, 


,350±0, 


.075 


0, 


.402±0. 


,051 


0, 


,309±0 


.109 


0. 


,205±0, 


,030 


0, 


,244±0, 


.016 


0, 


.223±0. 


,010 


OIII2 





.525±0, 


,114 


0, 


,604±0. 


,296 


0, 


.508±0, 


,300 


0, 


,384±0, 


.101 


0, 


.281±0. 


,086 


0, 


,223±0, 


.218 


0. 


,226±0, 


,128 


0, 


,278±0, 


.100 


0, 


.141±0. 


,043 


Fe5015 


1, 


.606±0, 


,361 


1, 


,668±0, 


,327 


1, 


.446±0, 


,310 


0, 


,789±0 


.127 





.767±0, 


,070 


0, 


,814±0 


.105 


0. 


,520±0, 


111 


0, 


,483±0 


.009 





.544±0, 


,066 


Mgi 


0, 


.021±0, 


,012 


0, 


,020±0, 


,003 


0, 


,020±0, 


,005 


0, 


,009±0 


,003 


0, 


,010±0, 


,001 


0, 


,009±0, 


,003 


0, 


,006±0, 


,001 


0, 


,006±0 


,000 





,007±0, 


,001 


Mg2 


0, 


.022±0, 


,006 


0, 


,022±0, 


,006 





,024±0, 


,006 


0, 


,012±0 


,002 


0, 


,015±0, 


,004 


0, 


,010±0, 


,002 


0, 


,008±0, 


,001 


0, 


,007±0 


,001 





,007±0, 


,001 


Mgb 


0, 


.742±0, 


,265 


0, 


,816±0, 


,186 


0, 


,779±0, 


,048 


0, 


,377±0, 


,027 


0, 


,406±0, 


,137 


0, 


,391±0, 


,067 


0, 


,237±0, 


,013 


0, 


,243±0, 


,031 


0, 


,242±0, 


,015 


Fe5270 


0, 


.931±0, 


,204 


0, 


,904±0, 


,367 


0, 


,910±0, 


,141 


0, 


,405±0, 


,072 


0, 


,447±0, 


,072 


0, 


,480±0, 


,125 


0, 


,260±0, 


,050 


0, 


,275±0, 


,014 


0, 


,304±0, 


,018 


Fe5270s 


0, 


.634±0, 


,154 


0, 


,644±0, 


,291 


0, 


,665±0, 


,082 


0, 


,277±0, 


,048 


0, 


.277±0, 


,043 


0, 


,318±0, 


,124 


0, 


181±0, 


,038 


0, 


,294±0, 


,158 


0, 


,204±0, 


,010 


Fc5335 


0, 


.882±0, 


,288 


0, 


,905±0, 


149 


0, 


,817d::0. 


,208 


0, 


,506±0 


,068 


0, 


.461±0, 


,048 


0, 


,390±0, 


,050 


0, 


,284±0, 


,050 


0, 


,406±0, 


,080 


0, 


,318±0, 


,038 


Fc5406 


0, 


.669±0, 


,252 


0, 


,664±0, 


,299 


0, 


,734±0, 


,094 


0, 


,339±0 


,072 


0, 


.376±0, 


,075 


0, 


,284±0, 


,078 


0, 


191±0, 


,023 


0, 


,243±0, 


,004 


0, 


,267±0, 


,089 


Fe5709 


0, 


.553±0, 


,152 


0, 


,484±0, 


144 


0, 


,539±0, 


,109 


0, 


,237±0 


,013 


0, 


,234±0, 


,037 


0, 


,216±0, 


,035 


0, 


154±0, 


,030 


0, 


,167±0, 


,021 


0, 


,183±0, 


,018 


Fc5782 


0, 


.476±0, 


,102 


0, 


,395±0, 


,071 


0, 


,463±0, 


,087 


0, 


,196±0, 


,080 


0, 


,235±0, 


,030 


0, 


,175±0, 


,035 


0, 


139±0, 


,006 


0, 


,146±0, 


,010 


0, 


,131±0, 


,023 


NaD 


0, 


.601±0, 


,172 


0, 


,530±0, 


141 


0, 


,644±0, 


,112 


0, 


,271±0, 


,037 


0, 


,353±0, 


,083 


0, 


,260±0, 


,037 


0, 


175±0, 


,028 


0, 


,225±0, 


,013 


0, 


,191±0, 


,032 


TiOi 


0, 


.015±0, 


,004 


0, 


,015±0, 


,004 


0, 


.019±0, 


,005 


0, 


,008±0, 


.001 


0, 


.008±0, 


,001 


0, 


,007±0, 


.001 


0, 


,005±0, 


,001 


0, 


,005±0, 


.000 


0, 


.005±0, 


,001 


Ti02 





.013±0, 


,001 


0, 


,013±0. 


,002 


0, 


.015±0, 


,003 


0, 


,006±0, 


.001 


0, 


.007±0, 


,001 


0, 


,006±0, 


.001 


0. 


,004±0, 


,001 


0, 


,004±0, 


.000 





.005±0, 


,000 



^mcan S/N of SDSS g, r and i-band 

'^thc best fitting amplitude A of the considered line and the level of 
scatter rN in the residuals of the gandalf fit on Ha 
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weak Ha intermediate Ha strong Ha 




4000 4500 5000 5500 6000 6500 7000 10 20 30 40 50 60 10 20 30 40 50 60 10 20 30 40 50 60 

Wavelenglh[3l] S/N S/N S/N 

Fig. 1. — Error estimation for absorption-line indices. The observed spectrum of a sample galaxy is shown in the left three panels, 
together with its best gandalf fit (red lines). In each of these panels, the inset shows the SDSS spectrum and our fit in a specific index 
region (for H/3, Fe5015 and Mgb from top to bottom, respectively, and as indicated by the dashed vertical lines in each main panel), as well 
as the flux density range of 100 synthetic realizations of this particular spectrum (grey bands) based on the formal SDSS uncertainties on 
the flux density and our gandalf fit. The right panels show the range (from top to bottom for the same indices as the left) when considering 
a sample of 405 objects (grey points) with varying quality of spectra, as traced by S/sN, and different strengths of nebular emission, as 
traced by the Ha line. From left to right, the red and green points highlight objects with increasing strength of Ha emission, with green 
points showing objects within specific S/sN ranges, where the median value of the errors on the absorption-line indices and the standard 
deviations are shown (crosses and vertical lines). 
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Fig. 2. — Quality assessment for continuum with examples inducing bad fits. This figure shows a quality assessment for a continuum 
region using randomly selected ~46000 objects. Red crosses on the central panel denote the median (lower) and 1 a (upper) distribution 
for each S/sN bin and orange solid lines fit (note that there is a demarcation line for 1 ct in the direction of larger rN/sN.). The vertical 
dashed lines indicate a specific bin and its Gaussian distribution has been inserted on the top side as an example. The color filled dots 
correspond to the colored left and right panels and clarify the trend for quality assessments. The black solid lines represent the observed 
spectra and the coloured lines are the fits. Furthermore, the top and bottom examples with gray fits are the major reasons for the bad fits. 
These telluric contaminated spectra (top) and Broad Line Regions (bottom) are marked with black crosses on the central panel. A minor 
reason which brought about bad fits is also denoted by red and green filled dots. 



formal errors. Such parameters are affected in different 
ways by data mismatches. For the steUar kinematics, 
and in particular for the velocity dispersion, the most 
crucial match is the stellar continuum across the entire 
wavelength range, whereas the accuracy of the emission- 
line parameters (and of the absorption-line indices that 
are affected by nebular emission) are more sensitive to 
the quality of the fit in more localized spectral regions 
around the emission lines. 

The quality of the data, as routinely estimated by di- 
viding the average level of the flux density (hereafter S, 
for signal) by the level of the formal uncertainties for 
the latter (hereafter sN, for statistical noise), does not 
guarantee a good fit. Artifacts introduced during data 
reduction may not be picked up by such an S/sN ratio, 
whereas a severe template mismatch to the stellar con- 
tinuum or the presence of additional components in the 
nebular spectrum (e.g., a Broad Line Region or Wolf- 
Rayet features that are not included in our standard fit- 
ting procedure) may lead to biased results, even in a 



case of excellent data. A more direct way to assess the 
quality of our model would be to compare the level of 
fluctuations in the fit residuals (hereafter rN, for resid- 
ual noise) to the expected statistical fluctuations, sN. A 
rN/sN ratio close to unity indicates a good fit, as this 
ration corresponds to a reduced x^, which is also close 
to 1. 

3.1. Continuum fit 

In order to assess the quality of the fit to the stellar 
continuum and help in isolating biased measurements for 
the stellar kinematics, we measured the rN/sN ratio in 
four different wavelength regions of equal width that were 
devoid of nebular emission. Specifically, we adopted the 
4500 - 4700A, 5400 - 5600A, 6000 - 6200A and 6800 
- 7OOOA wavelength intervals and combined the corre- 
sponding values for the rN/sN ratio into a single average, 
excluding the largest value in order to avoid a spurious 
contribution of continuum bands, where the SDSS spec- 
tra happened to be broken. 
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Fig. 3. — Quality assessment process for [Olll] A5007 with three examples. Left : The quality assessment plane drawn by a central 
emission region using randomly selected ~46000 objects. The three examples which have different emission line widths are shown as Ex. 1,2 
and 3. Middle : Quality assessment plane given by a typical passband. The lower and upper red filled dots are the median and la at 
each S/sN bin. The black dashed and solid lines trace each point. Moreover, [Olll] A5007 emissions for these three examples are shown on 
left sub-panels in corresponding colors. The black line represents the observed spectrum and the colored one the fit. Once we derived the 
median and Icr from the typical passband, we measured the Nct of every [O III] A5007 emission. The newly derived Icr arc also shown using 
the same colors. Right ; 1,000 objects shown in different colors depend on Ncr. The objects below the median arc indicated in gray. 



As shown in Fig. 2, for ^46000 randomly selected 
SDSS galaxy spectra, the distribution of the rN/sN aver- 
age values for the fit to the stellar continuum as a func- 
tion of the quality of the spectrum itself, as measured 
by the S/sN ratio, measured and averaged in the same 
continuum bands. The vast majority of the spectra were 
well matched by our model, with rN/sN close to unity. 
Our ability to match the data also appears to have im- 
proved as the quality of the spectra increased. In fact, 
for data of excellent quality (high S/sN), the level of 
fluctuation in the fit residual became even smaller than 
the formal uncertainties in the flux density, suggesting 
that the latter errors must under-predict the real level 
of statistical fluctuations in the SDSS spectra in such 
high-quality regimes. 

Fig. 2 also shows a few representative examples of 
gandalf fits to the the SDSS spectra. As stated above, 
the quality of the data themselves does not guarantee 
success or failure in the fit. It is possible to have formally 
good fits, both in cases of data with good and relatively 
poor values for the S/sN ratio (orange and blue exam- 
ples), whereas our model did not necessarily fail only in 
a case of poor data (red example), but also when the 
quality was very good (green example). 

In general there are three reasons that our objects show 
rN/sN ratios that are significantly above unity. Most 
often a poor fit to the continuum is the result of the 
presence of a conspicuous Broad line region (grey exam- 
ples, lower panels), which was not noticed by the SDSS 
pipeline in the process of separating the galaxy from the 
QSO spectra. Although the presence of broad line re- 
gions (BLR) can be accounted for by gandalf (see § 3.3), 
in the framework of our standard recipe that includes 



only narrow- line regions, the presence of an additional 
component leads to a considerable mismatch outside the 
wavelength regions near the Ha and H/3 lines. A second 
kind of failure occurs for objects in which the adopted 
stellar library appears to be inadequate in matching the 
observed SDSS spectrum (dubbed as "Bizarre" in Fig. 2 
and already shown in the red and green examples) . This 
occurs in a very small fraction of objects ; however, it 
often appears the spectra have been contaminated by 
features in the sky foreground. Finally, for a third kind 
of spectra, a high rN/sN value is the result of telluric at- 
mospheric features at the red end of the spectra. When 
such features are present, the SDSS formal errors under- 
estimate the fluctuations around the stellar continuum 
introduced by the telluric features, biasing the rN/sN 
value in the 6800-7000 continuum band to high values. 
We note however, that in this case the fit to the rest of 
the spectrum is not necessarily poor. 

Having ascertained that the vast majority of our mod- 
els showed good fits for the stellar continuum of the SDSS 
spectra, we could assign, to any object with rN/sN above 
unity, a probability of being a true outlier (in the distri- 
bution shown in Fig. 2), and hence a galaxy to which 
the stellar continuum has been poorly matched. For this 
we determined, as a function of the quality of the spec- 
trum, the median of the rN/sN distribution and corre- 
sponding outlier-resistant standard deviation (cr). We 
used the rN/sN distributions in S/sN bins containing 
at least 100 objects (see the inset of the central panel 
of Fig. 2), which we determined using the Voronoi bin- 
ning scheme(Cappellari & Copin 2003), and fitted the 
resulting median and sigma values (red crosses) with a 
polynomial function (orange lines). 
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We then computed the distance of each point from the 
median Unc in Fig. 2, expressing its distance in units of 
sigma, Ncr. Since the rN/sN distributions in each S/sN 
bin were well matched by log-normal functions, such No- 
value readily gave the probabihty of a point being a true 
outher. For example, Ncr=3 indicated that there was a 
99.73% probabihty that the object was a true outher and 
only a 0.27% probability of being the result of random 
fluctuations. 

3.2. Nebular fit 

Similar to the case of the stehar continuum, we needed 
to check the quality of the fit with the emission lines in or- 
der determine possible biases in the physical parameters 
we wanted to measure, such as the line width or fluxes. 
Even when the nebular emission was clearly detected, 
it was possible that the underlying assumptions in our 
modeling approach were not met, due to the presence of 
a broad line region or to non-Gaussian line proflles, for 
example. 

In order to assess the quality of the emission-line flt 

using the previous approach, it was important to give as 
miich weight as possible to the deviations that arose only 
in the spectral regions around the; nebular emission, and 
to adopt, for the wavelength interval where the rN/sN 
ratio was computed, a width that scaled with the width 
of the lines. Indeed, if we were to evaluate the rN/sN ra- 
tio over a fixed number of sampling elements, at a given 
fit quality the deviations from the fits of lines of very 
different widths would contribute unevenly to the rN/sN 
budget, thus limiting our ability to spot bad fits in the 
case of narrow lines. For each emission line we thus com- 
puted rN/sN around the line within a passband of width 
equal to five times the gas velocity dispersion cTgas. To 
then assess the quality of an emission-line fit in a given 
object, we normalized the distance from the median line 
in the rN/sN vs. S/sN diagram using the UrN/sN value 
(at a given S/sN) that would be found if all of the rN/sN 
and S/sN ratios had been measured with passbands of 
identical widths corresponding to the value of the cxgas 
observed in the object of interest. Such a i7rN/sN was 
found by simply rescaling the (TiN/sN values obtained us- 
ing a fixed bandwidth, corresponding, for instance, to the 
typical value of Cgas found in the SDSS galaxies 90 
km s~^). For a good fit, the amplitude of the fluctua- 
tions around the median value for the rN/sN ratio (at a 
given data quality or S/sN value) scaled simply with the 
square root of the number of sampling elements used in 
computing the rN/sN ratio. 

There was a second complication in assessing the qual- 
ity of the emission-line fit. Even for systems with similar 
line widths, evaluating the amplitude of the statistical 
fluctuations in the rN/sN ratio at a given data quality 
was difficult due to the presence of a larger fraction of 
formally bad ffts compared to the case of the stellar con- 
tinuum. Put another way, it was more difficult to iso- 
late the tail of the outliers when estimating the standard 
deviation in the rN/sN histograms for the emission- line 
fit. To circumvent this problem, for every line that we 
wanted to check in terms of the quality of our fit, we de- 
fined a continuum passband in its immediate vicinity and 
used the rN/sN vs. S/sN distribution for this passband 
to determine the reference rN / sN median and sigma pro- 
files that were actually used to assess the quality of each 



fit. As previously mentioned, such profiles were derived 
by adopting for the continuum passband a single width, 
which was set to five times the typical cxgas value in the 
SDSS galaxies, or approximately 90km s~^. 

Fig. 3 illustrates how this procedure works for the spe- 
cific case of the [O iii] A5007 emission line. The left panel 
shows, for the same ~46000 galaxies shown in Fig. 2, the 
rN/sN vs. S/sN distribution for the [Oiii] A5007 fit, 
adopting for each object a bandwidth equal to five times 
CTjQjjjj . In this panel, we show the positions of three ob- 
jects with similar values of the rN/sN ratio around the 
[O III] A5007 line, but with different ctjQjjjj values, as in- 
dicated by the size of the symbols. 

If the median and Uj-n /sN profiles could be drawn in this 
panel and used to assess the quality of the [O ill] A5007 fit 
in these objects, we would have been led to believe that 
the goodness of the emission-line fit degrades as we move 
from left to right, to higher S/sN values. This is not quite 
the illustrated by the central panel. This panel 

shows the rN/sN vs. S/sN distribution obtained in the 
continuum passband adjacent to the [O III] A5007 line, 
centered at 5006. 7A and adopting a width equal to five 
times 90kni s^^. The dashed and solid black lines show 
the median rN/sN values and the level where the objects 
are one UrN/sN away from the median, respectively (us- 
ing the same method as for the continuum fit). For lines 
that are considerably broader than 90km s~'^ the statisti- 
cal fluctuations around the median ought to be consider- 
ably smaller since a broader passband containing a larger 
number of sampling elements was adopted. Conversely, 
the quality of the fit for much narrower lines should be 
compared against a higher bar. This is illustrated by the 
coloured solid lines in the central panel of Fig. 3, which 
shows the level where objects with the same CTjQjjjj of 
the three considered examples should lie when found one 
(Tj.]yf/gN away from the median rN/sN line. When such 
yardsticks are used to express the distance from the me- 
dian in terms of values, although the object with 
highest S/sN ratio (green point) lies further away from 
the median line, due to its narrower [Oiii] A5007 line, 
the quality of its gandalf fit was not much worse than in 
the case of the galaxy with the worse data (lowest S/sN, 
light blue) and was much better than for the object with 
the broadest [Olli] A5007 lines (light green). Using this 
quality-assessment scheme, the right-hand panel shows 
how the goodness of our emission-line fit can be ilhis- 
trated (albeit only for a small number of objects for the 
sake of clarity), by colour-coding each point in the rN/sN 
vs. S/sN diagram for the [Oiii] A5007 line according to 
the values derived in the rN/sN vs. S/sN diagram 
for the continuum band adjacent to [O ill] A5007. 

Typically, formally bad emission-line fits occur either 
in the presence of strong lines where even small devia- 
tions from a Gaussian profiles lead to high rN/sN val- 
ues, or in the presence of an additional BLR component. 
Non-Gaussian lines can have a doublc-horncd shape. For 
instance, in the case of the Balmer lines, tracing star- 
formation activity confined to a circumnuclear ring can 
show blue or red wings, presumably due to out or in- 
flows of gaseous material, or, can display the triangular 
Voigt profile that is typical of strong AGN activity. Our 
Gaussian fit could lead to biased measurements for the 
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Fig. 4. — Broad line region objects. The central panel shows how we picked out broad line regions from the entire DR7 galaxies. The 
median and Scr distributions are shown with red crosses and orange lines, which are shown on left and right panels. Left : The black and 
red lines represent the observed spectrum and fit, respectively. The two shaded regions on the top panel are the continuum and blue-side 
wing of the broadened emission region that was used to derive the flux ratio. The dashed horizontal bars represent the moan flux levels for 
the specifically defined continuum and wing. The inserted panel on the top is the fitting result, without the new prescription for the broad 
line region objects. Once we picked the broad line region objects from the entire DR7 galaxies, we re-fitted them to improve the fit. Right 
: The blue line is our fit. 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I 1 — I — I — I — I — I — I — I — I — I — n — I — I — I I .1 I — rr—T 




log([NII] X 6583/Ha) log([NII] X 6583/Ha) 

Fig. 5. — BPT diagnostic diagram including broad line regions. The objects morphologically classified into early type galaxies by FracDeV 
and the concentration index were drawn with filled contours. We determined the broad line regions from the "galaxies" in the SDSS DR7, 
which is shown in green dots in both panels. Many of the green dots moved to Seyfert regimes at the right panel, by applying a new 
prescription which masked the entire Balmer absorption series to a properly find continuum. Broad line regions, with the exception of 
widely broadened emission lines, need this new prescription to exactly fit their continuum. 
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position of the lines with bhie or red wings, and to an 
overestimation of the width of the true thermal compo- 
nent in lines with a Voigt profile. Our flux estimates, on 
the other hand, are not heavily affected by the presence 
of non-Gaussians profiles. This is not the case, however, 
for objects with an additional BLR component, in par- 
ticular when gandalf chooses to match the broacl lines 
rather than the narrow component of the Balmer spec- 
trum. Since the BLR fluxes can greatly exceed the flux of 
the narrow Balmer lines, when these objects are plotted 
on standard BPT diagrams, such as the [Nii] A6584/Ha 
vs. [O III] A5007/H/3 diagram, they can be easily spotted 
in regions with very low values for the [Nil] A6584/Ha 
or [Olll] A5007/H/3 ratio that are not occupied by other 
galaxies with well-matched nebular spectra from AGNs 
or star- forming regions (refer to the following section). 
Conversely, no object with a non-Gaussian profile ends 
up in similarly strange regions of the BPT diagrams. 

3.3. Broad Line Region 

In the previous two subsections, we argued that the 
presence of a broad-line region frequently causes a bad 
fit to both the emission lines and the stellar continuum, 
in particular near the Ha and H/3 lines. Yet, provided 
that the presence of a BLR is automatically spotted, the 
impact of such a component on the pPXF extraction of 
the stellar kinematics and the subsequent gandalf fit to 
both the stellar continuum and nebular lines can be easily 
dealt with. In fact, one needs only to add a second series 
of broader Balmer lines in the gandalf setup to place a 
broader mask around to Balmer lines during the pPXF fit 
and then include the BLR component in the gandalf fit. 

In order to recognize the presence of BLR, we looked 
for a "shoulder" in the continuum region adjacent to the 
Ha-|-[Nii] A6584 region, which would otherwise be rela- 
tively flat. For this, we computed the ratio between the 
flux levels observed in two spectral regions of equal width 
on the blue side of the Ha+[Nll] A6584 blend. The first 
was placed as close as possible to the [N ii] A6548 (avoid- 
ing any fiux from this line even for the broadest narrow- 
line region, ~ 300km s~^), and the second was located 
just beyond the extent of the broadest BLR (centering 
it on the blue side of the 6594A absorption line that is 
predominantly due to Calcium and Iron). Then, similar 
to our quality-assessment procedures, once such flux ra- 
tio was measured in a given object, we simply checked 
to determine if it was exceptionally high compared to 
the distribution of other flux ratios in spectra of simi- 
lar quality. Fig. 4 shows two examples of spectra with 
BLR components and how our procedure isolated them. 
This diagram juxtaposes the quality of the spectra in the 
Ha-h[Nii] A6584 region (using the S/sN ratio in the 6000 

62GGA window, § 4.1) and the flux ratio defined above. 
This figure also shows the gandalf fit according to our 
standard prescription and, after adding a BLR compo- 
nent, and illustrates how our standard gandalf fit docs 
not always fail while choosing to match the BLR region 
(as already noticed in § 4.2). 

To assess the efficiency of our criterion in spotting BLR 
components, we visually inspected the spectra and the 
residuals of the standard gandalf fit for 10,000 randomly 
chosen objects, labeling those where we recognized the 
presence of a BLR. Adopting a 3a cut to isolate excep- 



tionally high values for the previously defined flux ratio, 
our procedure selected 83% of the BLR that we identi- 
fied by eye, while wrongly deeming likely the presence of 
a BLR in 0.8% of the objects where we had found none. 
Although these figures are quite encouraging, it must be 
acknowledged that only the strongest BLR regions can be 
identified if the S/sN is small(i.e., poor-quality spectra). 

Even then, such BLR components appear in a non- 
negligible fraction of the SDSS spectra of nearby galaxies, 
as overall we found 8,610 objects with a BLR component 
out or the entire DR7 sample, or almost 1.3% of the 
galaxies. 

As mentioned in § 4.2, the impact of a failed gandalf 
fit due to a BLR can be crucial when classifying the 
nebular emission with standard BPT diagrams. In fact, 
when the gandalf algorithm opts to fit the BLR region 
within our standard prescription, the object will move 
to a strange region in the BPT diagram, below and to 
the left of the area occupied by the star-bursting systems. 
Adding a BLR component in the gandalf fit rectifies this 
problem, and usually moves the narrow-line components 
of these galaxies back to the regions of the BPT diagram 
that are occupied by active galactic nuclei, either Seyfert 
nuclei or low-ionization nuclear emission regions (LIN- 
ERs). Such a correction was particularly evident when 
considering a sample of objects that rarely showed signs 
of central star formation, such as early-type galaxies. 
Fig. 5 shows an area occupied by the early-type galaxies 
in our DR7 sample (as morphologically classified in § 7.2) 
in the standard [Nil] A6584/Ha vs. [Oiii] A5007/H^ di- 
agnostic diagram, as well as the position of those objects 
where we spotted the presence of BLR, before and af- 
ter adding a BLR component in the gandalf fit, in the 
left and right panels, respectively. For the vast major- 
ity of the objects with a BLR, our procedure failed to 
match the narrow IIa-|-[Nii] A6584 lines in the standard 
gandalf fit and they resided in the lower left side of 
the BPT diagram. If we properly dealt with the BLR 
component, however, they generally moved to the region 
occupied by Seyfert nuclei. 

4. DATABASE 

Our measurements of the stellar and nebular kine- 
matics, for the strength of both the emission lines and 
the stellar absorption features, for the diffuse and neb- 
ular reddening components, as well as our assessments 
of the quality of the fit to the stellar continuum and 
each of the emission lines, formed a rather complicated 
database. To facilitate distribution, we decided to orga- 
nize our pPXF, gandalf and absorption-line index mea- 
surements in IDL structures, which allowed us to group 
all the emission- and absorption-line quantities. Such 
structures were then saved in single FITS files for each 
of the objects in our public catalogue. These files can be 
queried for single or multiple objects by their SDSS ID, 
galaxy coordinate and redshift via a web interface, or can 
be downloaded in single archive. For the case of single 
objects, we also provide a second FITS file showing the 
gandalf fit to the object of interest. 

The output parameters in the IDL structures are orga- 
nized following the logical order of our fitting procedure. 
First we found the pPXF outputs, the stellar redshift z^, 
and velocity dispersion cr* (in km s"^). These are fol- 
lowed by the nominal quality of the spectra and of our flt 
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Fig. 6. — Comparison between our values for the central stellar velocity dispersion and those computed by the SDSS pipeline (upper 
panels) and those listed in the MPA-JHU catalogue (lower panels). Starting from an original sample of randomly drawn 46,000 objects 
from our catalogue, the left panels shows the comparisons with matches in the SDSS and MPA-JHU databases, whereas the central and 
left panels shows objects with no nebular emission (A/N< 1 for [Olll] A5007, Ha and [Nil] A6584) or with strong lines (A/N> 4 for the 
same lines, at the same time), respectively. In each panel the dashed lines show the one-to-one line and the solid lines show the la contour 
for the plotted data. 



in the continuum, as quantified by the S/sN and rN/sN 
ratios introduced in §4.1, which relate to the probabihty 
that the fit is biased by some artificial or unaccounted 
feature, as given by the parameter. Negative N^. 
values occur when the rN/sN ratio are below the me- 
dian value, and always indicate a good fit. Next are the 
gandalf emission-line measurements, which we provide 
for all of the lines listed in Tab.l, in addition to the H/3 
line. For each line we list the observed and de-reddened 
fiuxes (in 10~"'^''ergs~^ cm~^ A~^), the equivalent width 
of the line (in negative A values), the A/N ratio that 
relates to how well a fine is detected (§ 2.2), and the N^- 
value (as defined in a §4.2). The strength of the lines 
and our ability to match them is then followed by the 
redshift and width of the recombination and forbidden 
lines (labelled as zhq, crna, ■Z[oiii]j '''[oiii]) well as by 
the value of the interstellar reddening Ay, which includes 
a second, nebular component Ay^neb, when the Balmer 
decrement is detected (§ 2.2). The presence of a BLR 
is also fiagged at this point. Although no detection cut 
based on the A/N values was applied to the emission- line 
measurements, these should be deemed unreliable when 
little or no emission is detected. Last, but not least, 
are the absorption-line measurements for each of the in- 
dices listed in Tab. 2, measured in A. It is most likely 
that only the objects with SDSS of the best quality (e.g., 
S/sN > 30) would prove useful. 
All the previously listed redshift, velocity dispersion, 



fiux, equivalent width and reddening values are also ac- 
companied by their corresponding formal uncertainties. 
In some cases, due to the presence of sky lines or gaps 
and artifacts in the SDSS spectra, the values of the pa- 
rameters can diverge, in which case they are flagged. 

We provide our entire database (line 
strengths, quality assessments, fitting SED) at 
http : / / gem . yonsei . ac . kr/ ossy/. 

5. COMPARISONS WITH OTHER PUBLIC 
DATABASE 

There are two existing databases of SDSS DR7 data 
with which we can compare our own measurements. The 
first consists of the SDSS DR7 pipeline outputs, and the 
second comes from the MPA-JHU release of DR7 spec- 
tral measurements^. As in previous releases, the SDSS 
DR7 pipeline does not remove the nebular emission prior 
to the absorption-line measurements, so we expected bi- 
ases in cases of considerable ionised-gas emission. On 
the other hand, by following the method of Tremonti et 
al. (2004), the MPA-JHU measurements were based on 
a procedure that is similar to ours, whereby the nebu- 
lar contribution to the spectra and the corresponding 
absorption-line infill was carefully treated. Yet, from 
the description of the MPA-JHU database, we deter- 
mined some technical differences between their proce- 
dure and the method adopted here which proved impor- 

^ Available at http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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Fig. 7. — Comparison between our values for the emission-line width and flux and those from the SDSS pipeline, starting from the 
same 46,000 randomly drawn objects used for Fig. 6. From top to bottom we compared the measurements for the [O III] A5007, Ha and 
[Nil] A6584 lines, for the line widths (in the left panels) and the line fluxes (in the right ones). The left and right panels are further split 
according to the global strength of the emission (weak : 2<A/N<5 for [Olll] A5007, Hq and [Nil] A6584, strong : A/N>10 for the same 
lines, at the same time). In each panel, the dashed line shows the one-to-one line and the solid contour shows the Icr contour level for the 
plotted data. 



tant for interpreting possible discrepancies with our mea- 
surements. For instance, we measured the emission-hne 
fluxes while simultaneously matching the stellar contin- 
uum and the nebular emission with gandalf , whereas 
the IVIPA-JHU emission-line measurements were carried 
out on the residuals of a previous fit to the stellar con- 
tinuum, presumably while masking the regions affected 
by emission, which we did when extracting the stellar 
kinematics with pPXF. 

5.1. Stellar Kinematics 

With the data in Fig. 6, we began by comparing, for 
the same subsample of randomly drawn ^^46000 objects 
used in Fig. 2 and Fig. 3, our pPXF velocity dispersion cr* 
measurements with the corresponding values from both 
the SDSS DR7 and MPA-JHU databases. The overall 
agreement between the various measurements was good, 
although our cr* values appeared to follow more closely 
the SDSS DR7 data than the MPA-JHU. This occurred, 
in particular, at the Iow-ct^, regime, which is inhabited 
mostly by relatively faint spiral galaxies showing consid- 



erable amounts of nebular emission. In fact, at the very 
low (7^ end of our sample the presence of emission caused 
most of the SDSS DR7 measurements to fail, and whereas 
the IVIPA-JHU catalogue reports ct^, values for these ob- 
jects, in most cases their values exceed ours. The large 
scatter in this regime can be explained by considering 
the low S/sN ratio of the spectra for such faint spirals, 
but such a bias cannot be easily explained. It is wor- 
thy of note, however, that based on simulations similar 
to those presented in § 2.5, and featuring SDSS spectra 
for both very old and young stellar populations, we found 
the pPXF can recover ct* values that scatter evenly around 
the input values down to a^, = 50km s~^ and S/sN = 10. 

5.2. Nebular Emission Fit 

Continuing with the nebular fit, Fig. 7 and Fig. 8 
show how our values for the width and the flux of strong 
lines such as [O ill] A5007, Ha and [N ii] A6584 compare 
against the same quantities from the DR7 pipeline and 
the IVIPA-JHU catalog, respectively. When comparing 
our values to the DR7 and JVIPA-JHU measurements, we 
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Same as Fig. 7, but showing a comparison between the MPA-JHU values and those reported in this work. 



used observed and intrinsic line- widths, respectively, and 
in both cases the flux values that were still affected by 
reddening due to dust in the SDSS objects. We further 
corrected the DR7 fluxes for foreground galactic extinc- 
tion to bring them in line with our and the MPA-JHU 
measurements, and brought the MPA-JHU fluxes back 
to their original values before the renormalisation for ex- 
tended sources. 

Our line-width and flux measurements agreed fairly 
well with the DR7 measurements (Fig. 7), in particular 
when the lines were strong and thus the manner by which 
the stellar continuum was accounted for should have had 
only a very limited impact on the emission-line width and 
flux estimation based on the Gaussian fits. Still, even at 
these regimes our Ha flux values appear to exceed the 
DR7 measurements by a few percent. Furthermore, as 
we consider weaker lines, the way the Ha flux measure- 
ments compared to each other appeared to be systemati- 
cally biased, rather than being simply affected by a larger 
scatter due to larger relative errors on the measurements 
of weak lines in both the DR7 and our dataset, as seems 
to be the case for [Oiii] A5007 and [Nii] A6584. Our flux 
values became systematically larger than the DR7 mea- 
surements as the Ha lines become weaker, which is to 



be expected given that we accounted for the presence of 
the underlying stellar spectrum and that in this regime 
the strength of the Ha emission became comparable to 
its corresponding stellar absorption feature. 

No such bias as a function of the emission-line strength 
was observed when we compared our Ha flux values with 
the MPA-JHU measurements, which do account for the 
stellar continuum (Fig. 8), although our intrinsic Ha line- 
width values were larger than the MPA-JHU estimates. 
From the limited information provided with the MPA- 
JHU on-line catalog, it is difficult to understand the 
cause of such small deviations; however, in the strong- 
emission regime where the flux estimates should be con- 
sistent between all three databases, our values were in 
agreement with the DR7 values, whereas the MPA-JHU 
measurements were not. 

5.3. Absorption Line Measurements 

We conclude this section by presenting Fig. 9 and 
Fig. 10, where we compare our line-strength measure- 
ments with the values for the same Lick indices provided 
by both the DR7 pipeline and the MPA-JHU catalogue. 
The DR7 release webpages specify that the DR7 mea- 
surements were made on the original SDSS spectra, with- 
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Fig. 9. — Comparison between our values of the H/3, Fe5015 and Mgb indices (top, middle, lower panels) and those from the SDSS 
pipeline, starting from the same 46,000 randomly drawn objects used for Fig. 6. As in Fig. 6, the left panels show all of the objects present 
in both databases, whereas the central and right panels show only the objects with weak nebular emission or with strong emission-lines. 
For the H/3 and Fe5015 indices, we defined weak objects as those with A/N values below 1 for the Ha and [Olll] A5007 lines, and active 
objects as those with A/N> 10 for those lines. In the case of Mgb, we required intermediate A/N values for Ha, [Olll] A5007 and the 
[Nl] AA5198, 5200 lines simultaneously to pick objects with weak emission. Conversely, we required strong emission in all three lines in 
order to select objects where the [Nl] AA5198, 5200 contamination to the Mgb index was important. In each panel the dashed line shows 
the one-to-one line and the solid contour shows the la contour level for the plotted data, as in Fig. 6, in this figure the dashed lines also 
show the contours plotted in the previous panels (from left to right). 



out degrading their resolution to match the resolution of 
the Lick/IDS system. Furthermore, a first inspection of 
the objects with good quality spectra and that were com- 
pletely devoid of emission revealed that this was also the 
case for the MPA-JHU measurements and that in both 
cases no attempt was made to correct the line-strength 
measurements for the impact of kinematic broadening. 
The central panels of Fig. 9 and Fig. 10 show how, for 
passive objects, our uncorrected measurements were well- 
matched with both catalogues. 

When objects with emission are considered, however, 
the DR7 values (Fig. 9) were heavily biased for the 
Lick indices where an emission line could fall either in 
the index passband or within the adjacent regions that 
were used to estimate the level of the pseudo-continuum 
above the absorption line region. In fact, during the 
DR7 pipeline analysis, the nebular emission was not sub- 
tracted prior to the line-strength measurements. When 



the emission-line fell within the index passband, as hap- 
pened for the H/3 and [Oiii] A5007 lines in the case 
of the H/3 and Fe5015 indices (top and central panels 
of Fig. 9 and Fig. 10), the strength of the absorption 
line was under-estimated. When the emission fell in 
one of the continuum passbands, as happened for the 
[Ni] AA5198,5200 doublet in the case of the Mgb index 
(lower panels), the line-strength was over-estimated, as 
it was measured against an artificially enhanced pseudo- 
continuum level (see also Emsellem, Goudfrooij, & Fcr- 
ruit 2003). On the other hand, the MPA-JHU line- 
strength measurements tended to agree with ours in the 
presence of nebular emission, as they were also corrected 
for line-infill (Fig. 10, top and middle panels for H/3 and 
Fe5015). Yet, the MPA-JHU catalogue does not include 
several of the emission lines that are present in our cat- 
alogue, such as the [Ni] AA5198, 5200 doublet, which ex- 
plains why disagreement between the Mgb measurements 
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Fig. 10. — Same as Fig. 9 but showing a comparison between the MPA-JHU values and those reported in this work. 



still persists (Fig. 10 (lower panels)). 

As a final remark, it should be noted that the nec- 
essary correction for the impact of kinematic broaden- 
ing can be quite severe in the most massive galaxies 
with correspondingly high values for the central stellar 
velocity dispersion (Kuntschner 2004). This occurs, in 
particular, for the case of indices, such as Fe5015 and 
Mgb, that were typically used to estimate the metal- 
licity of the degree of a-enhancement of stellar popula- 
tions, for which not accounting for kinematic broadening 
would lead to an underestimation of their values. When 
we considered that the MPA-JHU measurements for the 
Mgb index could have been over-estimated in the pres- 
ence of [Ni] AA5198,5200 emission, we note that using 
the MPA-JHU values could lead to over-estimated val- 
ues of the Mgb/Fe5015 ratio, which is a useful gauge of 
Qf-enhancement (Worthey, Faber, & Gonzalez 1992). 

5.3.1. The impact of the [Nl] AA5198, 5200 doublet on 
element abundance estimates 

To give an example of the potential extent of this prob- 
lem, we considered two different subsamples: one com- 
prising spiral galaxies displaying nebular emission, and 
a second including only quiescent early-type galaxies. 



These objects were morphologically selected as described 
in the following section and, after applying our standard 
A/N > A and A/N < 1 requirement for strong lines in 
order to pick active and quiescent objects (in addition 
to requiring the detection or not of [Nl] AA5198,5200 
emission), we assembled 311 and 579 objects for spirals 
and early-type galaxies, respectively. Adopting an in- 
dicative luminosity-weighted age of 2 and 10 Gyr for the 
stellar populations of these two classes of galaxies, we 
compared the position of ours and the MPA-JHU mea- 
surements for the Iron and Magnesium sensitive <Fe>^ 
and Mgb indices to the grids of Thomas et al. (2003) 
to determine stellar population models of varying metal- 
licities and overabundance of a elements, such as Mag- 
nesium(Fig. 11). For the case of passive objects (right 
panel), the 1-cr contours for the distribution of ours and 
the MPA-JHU measurements overlapped well and indi- 
cated an a enhancement with respect to the solar values 
of the [a/Fe] ratio. On the other hand, for the case of 
spirals (left panel), not accounting for the presence of 
the [Ni] AA5198, 5200 lines led to an overestimation of 

This index is defined in Gorgas, Efstathiou, & Aragon Sala- 
manca (1990) as < Fe >= (Fc5270 + Fe5335)/2 



18 



Oh et al. 



4 r 



2 Gyr:active late-type 



this work 



MPA-JHU 




■1.35 



2 3 4 

Mgb[A] 



4 r 



3 r 



A 

^ 2^ 



1 7 



ot 



10 Gyr:quiescent early- type 

[a/Fe]=0.0 
0.3 




0.5 
0.35 



0.05 
[Z/H]=-0.35 



-1.35 



12 3 4 5 6 

Mgb[/1] 



Fig. 11. — Effect of [Nl] AA5198,5200 emission on the Mgb absorption-line index and connected estimates of the a-element abundance. 
The left panel shows a comparison of the position of a sample of spiral galaxies with strong nebular emission (shown by the Icr contour 
level of the data) in the Mgb vs. <Fe> diagram with the model grids of Thomas et al. (2003). The solid and dashed contours indicate the 
distribution of the data from our catalogue and from that of the MPA-JHU database. Similarly, the right panel shows the same quantities, 
but for a sample of quiescent early-type galaxies (see § 5.3.1 for more detail). The model grids assume a luminosity-weighted mean age of 
2 Gyr and 10 Gyr for the spirals and early-type galaxies, respectively. 



the Mgb/<Fe> ratio, and therefore an abundance of a 
elements in the central regions of spiral galaxies which, 
from our database, would seem to have solar abundances. 

6. GALAXY SPECTRAL TEMPLATES 

Large-scale spectroscopic surveys, such as SDSS, have 
made it possible to create spectral templates for galax- 
ies of various kinds and morphological classes, which can 
be used as inputs in SED fitting routines where optical 
template spectra are needed, such as in photometric red- 
shift codes (e.g. Stebbins & Whitford 1948; Baum, W. 
A. 1962; Bolzonella, Miralles & Pello 2000; Wolf et al. 
2004). Therefore, we present here template spectra of 
elliptical, SO, Sa, Sb, Sc, and edge-on late-type galaxies, 
in the wavelength range between 3,799A and 7,000A( see 
Figs. 12, 13 and 14). 

We began the construction of our stacked galaxy spec- 
tra by restricting ourselves to a volume-limited parent 
sample of objects with Af,. < -17.62 and in the narrow 
redshift range of 0.020 < z < 0.025. This ensured a 
sufficient quality for the spectra to enable analysis and 
stacking, and to ensure that galaxy evolution would not 
be a concern given than such a redshift slide covers only 
~ 67 Myr. 

We then applied a set of selection criteria based on 
SDSS estimates for the galaxy size, inclination and light 
concentration in order to obtain early- and late-type 
galaxies for which we could check the morphology by 
eye. To allow such a visual inspection, we restricted our- 
selves to sufficiently extended and face-on SDSS galax- 
ies, by imposing a lower limit on the angular size as 
isoAr > 40" and the axis ratio as isoBr/isoAr > 0.4. 
Furthermore, we initially split early-type galaxies from 
their late-type counterparts by looking for objects with 
concentrated light profiles that were best matched by a 
de Vaucouleur profile. For this we required values for the 



FracDeV parameter, which measures how well the light 
profile is matched by a de Vaucouleur law rather than by 
an exponential profile, above 0.95 (see also Schawinski 
et al. 2007) and values for the concentration index Cr 
{Cr — tqq I rso where rgg and rso are the radii containing 
90 and 50% of the light, respectively) above 2.6 (see also 
Shimasaku et al. 2001; Strateva et al. 2001). 

Through this pre-selection we obtained 2591 galaxies, 
which were all carefully inspected by eye and which in 
the end yielded 452 objects, subdivided in E, SO, Sa, Sb 
and Sc morphologies (as listed in Tab. 4). This included 
175 late- type and edge-on galaxies as well, which were 
selected in a similar way, but by imposing isoBr/isoAr < 
0.3. Our selection criteria is summarized in Tab. 4. 

With these robustly classified objects at hand, we pro- 
ceeded to combine their SDSS central spectra by stack- 
ing the best gandalf spectral fit for each of them, once 
each of our models was brought back to the rest frame 
and normalized in flux at 5,500A. In this process we also 
recorded, at each wavelength, the Icr scatter of the flux 
values contributed by each of the stacked spectra, in or- 
der to estimate how much the spectral energy distribu- 
tions of the galaxies entering our templates differed from 
each other. Fig. 12 presents all such spectral templates 
(except for the edge-on late-type spirals), which extend 
from 3,799A to 7,000A, and Figs 13 and 14 show each of 
the spectral templates by itself, together with SDSS im- 
ages for the twelve sample galaxies in each morphological 
class that we selected. 

As a concluding remark, we note that our stacked spec- 
tra encompass only the central regions of the objects we 
selected. Specifically, on average the 3"-wide SDSS fibers 
sample 20% of the light from the E and SO galaxies in our 
sample, and 5% for the late-type galaxies. For our E, SO, 
Sa, Sb, and Sc galaxies, the difference between their total 
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Fig. 12. — Stacked galaxy templates for the various morphological types listed in Tab. 4 (except edge-on late-type spirals), normalized 
at 5,500 A. 



TABLE 4 

The properties of the sample galaxies used for the SED stacking 



Morphology'' 




FracDeVg'^ 


FracDeVr"* 


FracDeVi'^ 




isoBr /isoAr^ 






Fractiorir^ 


Elliptical 


111 


l.OitO.O 


l.OitO.O 


l.OitO.O 


3.3it0.2 


0.8it0.1 


69.1itl9.9 


5.9ltl.8 


17.1it4.2 


SO 


85 


l.OitO.O 


l.OitO.O 


l.OitO.O 


3.3it0.1 


0.7it0.1 


68.4itl7.1 


5.6itl.5 


19.7it4.6 


edge-on late-type 


175 


O.OitO.O 


O.OitO.O 


O.OitO.O 


2.4it0.1 


0.2it0.0 


61.1itl8.7 


23.7it5.2 


8.7it2.6 


Sa 


28 


O.OibO.O 


O.OitO.O 


O.OitO.O 


2.0it0.1 


0.6it0.2 


57.6itl3.4 


25.2it4.9 


4.6itl.7 


Sb 


32 


O.OitO.O 


O.OitO.O 


O.OitO.O 


2.0it0.2 


0.6it0.2 


60.0itl6.9 


26.9it4.2 


4.9it2.6 


Sc 


21 


O.OitO.O 


O.OitO.O 


O.OitO.O 


1.9it0.2 


0.7it0.2 


50.7it8.2 


25.5it4.3 


4.5itl.2 



^ All galaxies wore visually inspected by eye 

b Number of stacked galaxy spectra from a complete volume-limited samples 

Mean FracDeV parameters with 1 cr standard deviation in g,r and i bands provided by SDSS pipeline. In cases of elliptical and SO, all 
parameters are greater than 0.95 while late-type galaxies have smaller than 0.05. 

Mean concentration index(C7. = ^90/^50) with 1 cr standard deviation in r band. The value 2.6 used here to classify early-(>2.6) or 
late- type (< 2. 6) galaxies. 

^ Mean ratio between isophotal major and minor axis with 1 cr standard deviation in r band. All classes have greater than 0.4 of axis 
ratio except for edge-on late-type(<0.3). 

^ Mean major axis size with 1 cr standard deviation in r band. All chosen samples are greater than 40 arcsec on its major axis size. 
' Mean effective radius with 1 cr standard deviation in r band. 

J Mean light fraction(per cent) with 1 cr standard deviation in r band. This parameter shows the light fraction within 1.5 arcsec fiber 
radius. 



magnitudes and their fiber magnitudes is 0.02, 0.01, 0.09, 
0.12 and 0.07 mag, respectively. Obviously, our tem- 
plates are more representative for early-types than for 
late-type. Even for late-types, the systematic trend along 
the morphological type is still visible in our templates, 
and the systematic offset from the total light properties 
is only of order 10%. 

7. SUMMARY 

We presented a new database of absorption and 
emission-line measurements for all galaxies within a red- 
shift of 0.2, for which there are central spectra in T^^ 
SDSS data release. We modeled the stellar and neb- 
ular components of such SDSS spectra, adopting pub- 
licly available IDL machineries, with the aim of extract- 



ing both reliable and reproduceable measurements for 
the stellar velocity dispersion, the strength of various 
absorption-line features, and for both the flux and width 
of the recombination and forbidden emission lines that 
were observed in the central region of our sample galax- 
ies. In an attempt to improve on existing databases of 
SDSS spectral measurements, we included in our fit to 
the stellar continuum both standard stellar population 
models and empirical templates. Our match for the neb- 
ular spectrum features an exhaustive list of both recom- 
bination and forbidden lines. Foreground dust Galac- 
tic extinction is also implicitly treated in our models 
as a simple dust screen component affecting the entire 
spectrum, whereas a second reddening component affects 
only the nebular spectrum. 
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Fig. 13. — Stacked templates and optical images for elliptical, lenticular and edge-on late- type galaxy morphologies. The left panel 
includes the Icr scatter estimated while stacking the spectra, with the aim of conveying how much difference there was in the shape of the 
spectral energy distribution of the galaxies used in each morphological template. As in Fig. 12, the templates were normalized at 5,500 A. 
Each image is 100.0" X 100.0" across, corresponding, on average, to a physical size of 100 kpc. 
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Our absorption-lino measurements compare well with 
the SDSS pipeline values in the absence of emission. Con- 
versely, the nebular flux and line-width measurements 
agree in the strong line regime. This is consistent with 
our superior separation of the nebular and stellar spec- 
tra. On the other hand, since The MPA- JHU catalogue is 
based on a similarly careful procedure, our measurements 
also agree when considering weak emission, although 
some emission lines were omitted in the MPA- JHU cat- 
alog. This is the case for the [Ni] AA5198,5200 lines 
in particular, which affect estimates for the a-element 
abundance of galaxies. 

Most importantly, our absorption- and emission-line 
measurements come with a quality assessment for the 
spectral fit that yielded them, in both the stellar con- 
tinuum and the nebular emission and across different 
wavelength regions of the SDSS spectra. Such a quality 
assessment also allows the identification of objects with 
either problematic data or with peculiar features that are 
not accounted for by the rather standard set of assump- 
tions that enter into our models. For instance, based on 
the quality assessment around the Ha and [Nil] A6584 
lines, we report here that approximately 1% of the SDSS 
spectra classified as "galaxies" by the SDSS pipeline do 
in fact require additional broad lines to be matched. 

We conclude our work by providing new spectral tem- 
plates for galaxies of different Hubble types, obtained by 
combining the results of our spectral fits for a subsample 
of 452 morphologically selected objects. As they repre- 
sent the average spectral energy distribution of galaxies 
of a given Hubble type, these templates should be well- 
suited for photometric redshift estimates based on broad- 
band fluxes as long as the target galaxies appear to have 
a relaxed morphology. The spectral variations with nar- 
rower wavelength ranges, due for instance to the impact 



of a different mass on the strength and width of stellar 
and nebular features, will be the subject of future inves- 
tigations based on larger parent sample for the stacked 
templates that could be selected, for instance, with less 
stringent morphological or redshift criteria. 

Finally, with the hope of further fostering the discovery 
potential of our database, we also made all of our spectral 
fits available to the community. 
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